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The nuclear spins of ionized donors silicon have become an interesting quantum resource due 
to their very long coherence times. Their perfect isolation, however, comes at a price, since the 
absence of the donor electron makes the nuclear spin difficult to control. We demonstrate that 
the quadrupolar interaction allows to effectively tune the nuclear magnetic resonance of ionized 
arsenic donors in silicon via strain and determine the two nonzero elements of the S-tensor linking 
strain and electric field gradients in this material to S n = 1.5 x 10 22 V/m 2 and S 44 = 6 x 10 22 
V/m 2 . We find a stronger benefit of dynamical decoupling on the coherence properties of transitions 
subject to first-order quadrupole shifts than on those subject to only second-order shifts and discuss 
applications of quadrupole physics including mechanical driving of magnetic resonance, cooling of 
mechanical resonators and strain-mediated spin coupling. 


The potential of computing based on quantum me¬ 
chanics has triggered a quest for a scalable quantum bit 
or qubit technology with local control and long coher¬ 
ence times. Recent studies on the nuclear spin of ionized 
phosphorus donors in silicon have demonstrated ultra- 
long coherence times X 2 exceeding tens of minutes even 
at room temperature m- Together with the advantages 
of modern silicon semiconductor technology, these times 
make such nuclear spins a prime candidate for the imple¬ 
mentation of a quantum memory [3]. To enable storage 
of quantum information on longer timescales, the imple¬ 
mentation of error correction is vital [4], requiring local 
control of single qubits, as well as two- or many-qubit in¬ 
teractions Such a local control can be achieved by the 
combination of the global use of nuclear magnetic reso¬ 
nance (NMR) with a local tuning by gates [6], which has 
proven challenging even for the neutral phosphorus donor 
system [7j. For ionized donors, the absence of the corre¬ 
sponding donor electron impedes any electronic control 
which would otherwise be mediated by the so-called hy- 
perfine interaction of electron and nuclear spins, leading 
to a near perfect isolation of the ionized donor nuclear 
spin. In contrast to phosphorus, however, the heavier 
group V donors with a nuclear spin I > 1/2 such as ar¬ 
senic possess a nonzero quadrupole moment which inter¬ 
acts with electric field gradients 0. Here, we show that 
this additional interaction can be manipulated by elastic 
strain applied to the host crystal and can very success¬ 
fully be used to mechanically tune the properties of ion¬ 
ized donor nuclear spins. This paves the way for scalable 
local addressing of these qubits, as strain can be applied 
locally by piezo-actuators on the scale of nanometers [9]. 
Similarly, mechanical resonators in the required MHz fre¬ 
quency range are routinely fabricated in silicon in the 
framework of microelectromechanical systems (MEMS) 


and should be able to realize spin-spin couplings medi¬ 
ated by phonons, ultimately generating two- or many- 
qubit gates H03- Furthermore, while the lack of a sta¬ 
ble silicon isotope with I > 1/2 has prohibited accoustic 
nuclear resonance studies into quadrupole interaction in 
this material, the experiments reported here finally al¬ 
low to measure the effectiveness with which electric field 
gradients are generated on substitutional lattice sites in 
silicon under the application of strain. 

The initialization of ionized donor qubits can be real¬ 
ized by processes of spin-selective ionization based on the 
hyperfine interaction of electron and nuclear spins in the 
neutral charge state, which in most cases allows for an 
electrical read out of the qubit at the same time [IIHI3]. 
In our experiments, we detect the resonance frequencies 
of the nuclear spins of ionized arsenic donors As + mak¬ 
ing use of the Pauli blockade formed in weakly coupled 
spin pairs. These pairs form between the neutral donors 
As 0 with an electron spin S = 1/2 and the paramagnetic 
oxygen-vacancy complex (SL1) in monocrystalline sili¬ 
con exhibiting S = 1 pams] [Fig. 1 (a)]. Because of the 
Pauli exclusion principle, the electron transfer from the 
donor level to the energetically lower SL1 level strongly 
depends on the spin symmetry of the two partners (par¬ 
allel and antiparallel). This leads to an effective polar¬ 
ization into the long-lived parallel configuration in the 
steady state. Changing the spin symmetry by manipu¬ 
lating the spin state of either partner promotes the charge 
transfer, which is detected as a resonant increase in the 
conductivity, attributed to an Auger-recombination-type 
process M- 

The samples used in this work are silicon with natu¬ 
ral isotope composition which was implanted with As + 
ions, resulting in a maximum As concentration of about 
7 x 10 19 cm -3 about 70 nm below the surface. The sam- 
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pies were not annealed in order to preserve the defect 
centers created upon implantation m, which leads to a 
low dopant activation ratio as confirmed by the observa¬ 
tion of hyperfine-resolved ESR peaks. The samples were 
contacted with an interdigit structure, biased with typi¬ 
cally 8 V and illuminated with a red (wavelength 635 nm) 
light-emitting diode. All measurements were performed 
at 8 K and at a fixed microwave frequency (9.74 GHz) 
in a dielectric resonator for pulsed ENDOR. After the 
application of mw pulses, a resonant transient increase 
of the conductivity of the sample is recorded with a fast 
digitizer card. 

The electron spin resonance (ESR) of As 0 splits into 
four different levels, in first order separated by the hy- 
perfine interaction of 198.35 MHz [TB]. Each ESR line 
corresponds to one of the four nuclear spin eigenstates 
mi = —3/2... 3/2 of the / = 3/2 75 As nucleus [Fig. 1 
(b)], allowing the selective ionization and read out of the 
nuclear spin state of the ensemble. In combination with 
radio frequency (rf) pulses, this allows us to perform elec¬ 
trically detected electron nuclear double resonance (ED 
ENDOR) mm, as shown schematically in Fig. 0(c). 
At the beginning of the pulse sequence, a selective mi¬ 
crowave (mw) pulse lifts the Pauli blockade and thus 
promotes the ionization of As 0 in a particular nuclear 
spin state, here m/ = +3/2 ( selective ionization). After 
the fast antiparallel recombination (on a timescale of 5 
fi s [15]), the ionized donor spins are highly polarized and 
can be manipulated by rf pulses ( NMR control ), changing 
the nuclear spin state. This leads to a nuclear polariza¬ 
tion after the slower parallel recombination (600 /xs) has 
ionized the remaining donors. After the whole ensemble 
is transferred to the neutral charge state via illumination 
by a light-emitting diode (LED) and electron capture, 
the polarization is measured by application of an elec¬ 
tron spin echo sequence and detection of the ensuing cur¬ 
rent transient through the sample (electrically detected 
magnetic resonance, EDMR) (read out) [T9H2T] . When 
ionization and detection are performed on the same reso¬ 
nance (mi = 3/2), the observed EDMR signal is reduced 
in case of successful NMR manipulation. 

The spin states of an ionized arsenic donor As + in sili¬ 
con are subject to the nuclear Zeeman interaction, which 
splits the energy levels in an external magnetic field Bq 
(assumed to be oriented in the z direction), and the in¬ 
teraction between the quadrupole moment q and an elec¬ 
tric field gradient Vij = dEi/dj , where Ei are the elec¬ 
tric field components and i,j are the coordinates x,y,z. 
The total Hamiltonian describing the nuclear spin com¬ 
ponents I x , I y , and I z can be written as 

H n = -hf 0 I z + V vWi + W - ^ 4 ] ’ 

id 

where /o = /a n g n Bo/h is the Larmor frequency of the 
nuclear spin, /x n and g n are the nuclear magneton and the 
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FIG. 1. Manipulation and electrical detection of 75 As nuclear 
spins based on spin-dependent recombination, (a) Charge 
states As 0 and As + of an arsenic donor with nuclear and 
electron spin states denoted by green and blue arrows, re¬ 
spectively. The ionization is governed by a Pauli blockade 
of the donor electron spin with a paramagnetic defect (SL1). 
(b) Schematic spectrum showing the enhanced conductivity 
caused by electrons ejected to the conduction band during the 
ionization process nu. ( C ) Evolution of the charge and spin 
states of the 75 As donors during the ED ENDOR sequence. 
See text for details. 


nuclear g-factor, resp., h is Planck’s constant and e is the 
elementary charge [8]. For a cubic crystal symmetry as in 
unstrained Si, the electric field gradients cancel out and 
the spin states are separated by the Zeeman interaction 
only. In this case, the frequencies of the three allowed 
NMR transitions (Am/ = ±1) are degenerate [Fig. [2] 
(a)]. 

To measure the NMR frequencies of As + , we record ED 
ENDOR spectra at four magnetic fields, corresponding 
to the ESR lines of As 0 [Fig. 2 (c)], and thereby detect 
the relative occupancy of the four nuclear spin states of 
As + . When monitoring these occupancies while sweeping 
the frequency of the rf pulse, a reduction is expected on 
resonance with the NMR transitions. Therefore, after 
ionization of, e.g., the m/ = 3/2 state, the signature of 
the transition m/ = 3/2 «+ 1/2 is anticipated, while for 
ionization on m/ = 1/2 the measurement is sensitive to 
two transitions, m/ = 3/2 o 1/2 and m/ = 1/2 —1/2. 

Experimentally, in unstrained samples only one peak is 
observed in each of the four NMR spectra [Fig. 2 (d)]. 
This is in agreement with the expected degeneracy of the 
three NMR transitions for zero quadrupole interaction 
[Fig. 2 (a)]. From the observed resonance positions g n = 
0.9558(2) can be extracted, corresponding to a chemical 
shift of —0.40(2)% with respect to the value for free As 
nuclei [22] . 

A uniaxial strain e± applied to the cubic Si crystal 
reduces the cubic symmetry and results in an electric 
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FIG. 2. Demonstration of the strain-induced tuning of ion¬ 
ized donors As + . (a) Effects of the first- and second-order 
quadrupole interaction on the magnetic sublevels of As + . (b) 
Generation of strain via the differences in the thermal expan¬ 
sion coefficients of Si and a sapphire substrate. As donors 
are generated by shallow implantation, (c) EDMR spectrum 
of As 0 in Si. (d) NMR transitions of As + in unstrained Si 
after ionization on each of the four As 0 EDMR peaks, (e) 
Corresponding NMR transitions for As + in strained Si show¬ 
ing first-order quadrupole shifts. The inset shows a high- 
resolution spectrum revealing the second-order quadrupole 
shift of the mi = 1/2 —1/2 transition. 


field gradient V± which interacts with the arsenic nuclear 
spin, giving a nonzero quadrupole energy Ji/q = V±eq. 
The change in resonance frequency is described by (i) 
a first order shift A/gi oc /g of the “outer” transitions 
mi — 3/2 1/2 and —3/2 —1/2, and (2) a second or¬ 

der shift A/q2 oc /g//o of all four resonance lines, includ¬ 
ing and most easily observed on the “inner” transition 
mi = 1/2 —1/2 [Fig. 2 (a)]. In our proof-of-principle 

experiments, strain is applied by cementing a thin silicon 
sample to a sapphire substrate which results in uniax¬ 
ial strain e± at low temperature, caused by the different 
thermal expansion coefficients of the materials [Fig. 2 
(b)]. A precise in-situ measurement of the resulting strain 
is provided by the change of the hyperfine interaction of 
electron and nuclear spin of the neutral donor in the case 
of strain applied in the [100] direction of the crystal {23j. 
We find an expansion of e± ~ 2.9 x 10 4 . This corre¬ 
sponds to an in-plane compression of ey ^ —3.8 x 10 -4 , 
in agreement with the thermal expansion coefficients for 
Si and c-plane sapphire [24] [25]. Due to the high repro¬ 
ducibility of the strain observed in all of the (100) Si- 
sapphire stacks investigated and the isotropy of the ther¬ 
mal expansion in Si, we assume that the applied strain 
is equal in (111) stacks, to which this method to mea- 




FIG. 3. Angular dependences of the first-order (a) and 
second-order (b) quadrupole shifts on the relative orientation 
of the magnetic field Bo and the electric field gradient V± 
generated in [111] (black) and [100] (blue) direction. 


sure €j_ cannot be applied because of the symmetry of 
the conduction band minima [23] . Figure 2 (e) shows 
that for (111) stacks, the first order shift is clearly ob¬ 
served, separating the peaks corresponding to inner and 
outer transitions. The outer transition resonance peaks 
are shifted by well above one linewidth, while the sharper 
peaks of the inner transition are barely moved with re¬ 
spect to fo on the frequency scale of Fig. 2 (e). The 
lineshapes of the sharp lines close to fo in Fig. 2 (d) and 
(e) are dominated by a pattern reflecting the distribu¬ 
tion of frequencies given by the Fourier transform of the 
400 fi s NMR square pulse used, while the shapes of the 
resonances observed for mi = ±3/2 are asymmetrically 
broadened. For the nominally unstrained samples, this 
is most likely due to strain induced by the contact struc¬ 
ture m, the additional asymmetry of the broad lines in 
Fig. 2 (e) is attributed to an inhomogeneity of the strain 
generation. The weak resonances at / — /o = 0 in the 
top and bottom spectra most likely originate from a dy¬ 
namic crosstalk, where changes in the polarization of the 
mi = 1/2 and —1/2 states influence the polarization of 
the monitored states via relaxation. 

When using longer rf pulses, a higher spectral resolu¬ 
tion can be reached. For a pulse length of 6 ms, the sec¬ 
ond order shift to the inner transition becomes evident, 
as shown in the inset of Fig. [2] (e), where the resonance 
line is shown with (black trace) and without strain (red 
trace) in direct comparison. Since the linewidth of this 
transition is smaller by more than two orders of magni¬ 
tude compared to the outer transitions, the resonance is 
shifted by more than one linewidth as well, even though 
the shift A/g 2 is significantly smaller in absolute num¬ 
bers. Note that the quadratic dependence of A/q 2 on 
/g will reduce the negative influence of small fluctua¬ 
tions for this shift. To further confirm our interpretation 
of the observed shifts of the resonances as resulting from 
quadrupole interactions, the characteristic dependence of 
both first and second order shifts on the angle $ between 
strain axis and magnetic field was measured [Fig. 3 (a) 
and (b)]. When treated as a perturbation of the Zeeman 
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interaction, the expected quadrupolar shifts for / = 3/2 
are given by [ 8 ] 


-(3 cos 2 $ — 1 ) 


and 


A/qi = & 

A/q 2 = - cos 2 i9)(9cos 2 i? - 1 ) , 

64/o 


which fits the observed angular dependence very well 
[solid lines in Fig. 3 (a) and (b)]. 

The quadrupole shift is caused by an anisotropic elec¬ 
tric field gradient which interacts with the arsenic nu¬ 
clear spin and is proportional to the applied strain, 
Vi = Sij • €j (using Voigt notation). The tensor 
describing the proportionality in cubic crystals has two 
nontrivial components, S±\ and S 44 , which quantify the 
effect of strain (along the principle axes) and shear, re¬ 
spectively m ■ They can be investigated separately by 
measuring samples strained in [ 100 ] and [ 111 ] crystal di¬ 
rections. When fitting the angular dependence functions 
to the experimental data for first and second order, an 
/q of 76(5) kHz and 255(3) kHz can be extracted from 
the shifts in [100] and [111] samples, respectively. Hence, 
Sn = 1.5(3) x 10 22 V/m 2 and 644 = 6.8(20) x 10 22 V/m 2 
using q = 31.4 fm 2 for 75 As [ 22 ] . 

While there are no previous measurements of S in 
Si due to the absence of stable isotopes with nonzero 
quadrupole interaction, we note that the obtained val¬ 
ues are similar to those reported from acoustic nuclear 
resonance on 73 Ge in Ge [28], both with respect to the 
absolute values of Six and S 44 and to their relative size. 
This correspondence is expected due to the similar crys¬ 
tal structure of Si and Ge and the equal core electron 
configuration of As and Ge. Experiments on Sb and Bi 
will show how universal the ^-parameters obtained here 
are for silicon. 

To study the influence of the quadrupole interaction 
on the decoherence of the nuclear spins, we measure the 
coherence time X 2 using a spin echo sequence with an ad¬ 
ditional 7r/2 projection pulse which is inserted after the rf 
7 r-pulse [Fig. 4 (a)]. Due to the larger coupling strength of 
the first order effects, we expect a bigger influence on the 
coherence time observed on the outer transitions. The 
decay of the echo amplitude [black traces in Fig. 4 (b) and 
(c)] confirms the expected difference, giving X 2 ~ 44(5) 
ms for the inner and X \ ^ 3.7(3) ms for the outer tran¬ 
sition. These values are found to be mostly unchanged 
after the application of stress in our experiments. While 
X 2 for the inner transition is comparable to the coherence 
time of ionized phosphorus donors in silicon with natural 
isotope composition m, the coherence time of the outer 
transition is reduced, in accordance with the expected 
stronger coupling to the phonon bath. By application of 
a Carr-Purcell dynamical decoupling sequence [29] with 
32 refocusing pulses, X 2 is extended to X 2 = 275(10) ms 
and X 2 = 64(5) ms for the inner and outer transitions, 
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FIG. 4. (a) Pulse sequence used for coherence time measure¬ 
ments. Nuclear spin echo decays of the 3/2 ++ 1/2 (b) and 
1/2 —1/2 (c) subsystems of the As + nuclear spins. Colors 

denote different numbers of refocusing pulses n in the dynam¬ 
ical decoupling sequence as indicated in the insets, dashed 
and solid lines correspond to fits with superexponential de¬ 
cays exp[— (t/T 2 ) a ] with different exponents a. The insets 
show the extracted decay constants T 2 as a function of n. 


respectively. While the enhancement with the number 
of applied decoupling pulses n scales oc n 0,5 for the in¬ 
ner transition [inset of Fig. 4 (b)], as would be expected 
for 1/f noise [30] . a different scaling oc n 0,8 is observed 
for the outer transition [inset of Fig. 4 (c)], indicating 
that a different process dominates the decoherence in this 
case. This could be connected to, e.g., phonon number 
fluctuations j3J or field gradients introduced by charged 
defects. However, the stronger scaling with n suggests 
that, for a large number of pulses, this additional deco¬ 
herence mechanism can be made negligible. Even longer 
coherence times could most likely be reached using iso- 
topically purified silicon to overcome the interaction with 
29 Si nuclear spins in the host crystal [32] . 

Coming back to the application of the strain interac¬ 
tion to the tuning of qubits, with strains of 5 x 10 -5 
which can be realized by electrically addressable piezo¬ 
actuators [33] . shifts of the 3/2 1/2 transition of the 

order of 18 kHz could be realizable. Compared to the 
corresponding transition in our proof-of-principle sam¬ 
ples in Fig. 2 (e), this strain corresponds to a shift by 
one full linewidth. Going to smaller ensembles or sin¬ 
gle As + spins, strain inhomogeneities will be reduced or 
removed, allowing shifts of many linewidths and an ex¬ 
cellent selectivity, even without the use of isotopical en¬ 
gineering. Due to the strain connected to the oscillation 
of, e.g., a nano-mechanical beam, the quadrupole inter¬ 
action also renders possible the coupling of the nuclear 
spin to micro-mechanical modes. An estimate based on 
elasticity theory for a doubly clamped beam m with di- 
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mensions 1 x 0.1 x 0.1 fi m 3 results in a coupling constant 
g/27T « 18 Hz, which is of the same order of magnitude 
as the observed coherence time. Other possible applica¬ 
tion of strain interactions include the implementation of 
mechanical driving of nuclear spin resonance in Si and 
the nuclear spin-based cooling of mechanical oscillators, 
allowing to achieve higher displacement and force sensi¬ 
tivity [3H |35j, as well as the generation of squeezed spin 
states [361 . 
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